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Summary 

The steady discharge of a high-voltage corona is analyzed 
as a singular perturbation problem. The small parameter is 
the ratio of the length of the ionization region to the total gap 
length. By this method, current versus voltage characteristics 
can be calculated analytically. 


Introduction 

The equations governing the corona discharge which occurs 
when a thin electrode in a gas is raised to high voltage are 
nonlinear because the rate of ion production is a function of 
the electric field and the electric field is determined by the 
distribution of charge within the gas. However, a high electric 
field is required to ionize the gas, so ionization occurs only 
in a thin boundary layer around the cathode (Vatazhin and 
Grabovskii, 1983). It is natural to analyze corona discharge 
by a matched asymptotic analysis (Bender and Orszag, 1978). 
By doing so, we take advantage of the fact that the high electric 
field is generated by the small radius of the discharge electrode 
in order to find closed-form solutions for the corona. 

Our method of analysis will be applied to the negative corona 
occurring in a coaxial geometry although the method also 
works for positive corona. Many analyses of negative corona 
have treated pulse formation in a negative point— positive plane 
configuration. These consisted primarily of numerical 
simulations (Morrow, 1985; Davies and Donne, 1980). 
Indeed, the Trichel pulses produced at negative points are of 
considerable interest. However, the present analysis deals only 
with the more-or-less continuous discharge produced by a 
corona wire (Gallo, 1977). It is unclear at present whether 
the asymptotic method can be applied to the highly unsteady 
Trichel pulse corona. Research in that direction is currently 
in progress. 


Governing Equations 

The dependent variables in this problem are the space charge 
density and the electric field. Let L be a characteristic 
dimension of the overall geometry, a a typical radius of the 


cathode, and <p Q the anode potential. The magnitude of 
electric field required to produce ionization will be E h and 
the total current in the corona will be I 0 . Ultimately, l 0 is to 
be found. On a fundamental level, E ( would be an ionization 
potential divided by the mean free path. Then the following 
nondimensional variables can be formed: 
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where asterisks denote dimensional variables, q is the charge 
of an electron, n$ the number density of specie 0 ( + ions, 
- ions, or electrons), <p the electric potential, and a, and ct a 
are ionization and attachment coefficients (a* is Townsend’s 
first coefficient, which we take to be 0(1/6)). The b@ 
(/3 = + , - , or e) are the mobility coefficients appearing in 
Ohm’s law 


— ±b@E (2) 

where is the drift velocity and the minus sign is used for 
electrons and negative ions. The nondimensional b e equals 1. 
A small parameter equal to the scale of the ionization region 
divided by the scale of the overall geometry is given by 


= Vo_ 
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In order for discharge to occur, a must be 0(6). 

We will consider a negative corona with electrons generated 
at the cathode surfaces by secondary emission; electron 
generation by photoionization is assumed to be negligible. Ions 
are produced within the gas by electron impact. In 
dimensionless form, the steady-state equations for conservation 
of species and electric flux are (Vatazhin and Grabovskii, 
1983) 
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and e g is the gas dielectric constant. The value of the 
nondimensional parameter J determines I 0 when the other 
terms in equation (3b) are given. The electric field is 

E = - V<p (3c) 

In regions where \E\ is small compared to 1/e, a, will be 
zero. We will represent a, by the form 

a = j3e -1/{|£| (4) 


Because of the exponential dependence (eq. (4)), a, will be 
effectively zero, except within a distance of 0(e) around the 
cathode. Outside this region the production of ions by electron 
impact vanishes. Within the 0(e) boundary layer (or “sheath”), 
the spatial variation of <p is so rapid that the first of equations 
(3a) can be replaced by Laplace’s equation; the contribution 
from space charge is negligible compared to the geometrical 
divergence of the electric field. To develop the asymptotic 
analysis, we will consider the coaxial geometry shown in the 
following sketch: 



Asymptotic Analysis of Coaxial Discharge 


(See Meek and Craggs, p. 60; note that 0 and 1/e are 
proportional to the gas pressure.) Thus, where \E\ = 0(1), 
a, is exponentially small, and a a will be taken to be 
independent of E. 

Equations (3a) are to be solved with the boundary conditions 


<p = 0 on the cathode 

<p — l on the anode 


(5a) 


Because the reference current l 0 is arbitrary, it can be defined 
so that the net positive ion current is -1 on the cathode: 
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Inner Region 

Within the boundary layer we introduce the 0(1) scaled radial 
variable r = r/e. Then, with the notation np = mpVp, to 
lowest order, equations (3a) become 
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Cathode secondary emission occurs where positive ions knock in cylindrical coordinates. The solution to equation (6a) subject 

electrons out of the surface. If y is the coefficient of secondary to <p — 0 on r = a is 

emission, then 


yn + v + *h = —n e v e 9 n 


(5c) 



(7a) 


Also, 


n_ — 0 on the cathode 

n+ = 0 on the anode 


The constant A will be determined by matching with the outer 
solution. 

(5d) The solution to equations (6b) with x + = — 1 and 
= — yi r + on f = a is 
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£a,- dr 

ir e = 
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In equation (7b) a t is given by equation (4) with E = — Air 
from equations (3c) and (7a). Hence, 


Outer Region 

In the outer region, r = 0(1). Equations (3a) then become 


- — r— = J(n e + n_) 
rdr dr 
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On introducing equations (10), integrating equation (11), and 
matching to equation (7a), one finds 
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In order that 7r + in the inner region (eq. (7b)) match its 
value of 0 in the outer region, the relation 
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must be satisfied. This equation ultimately determines the value 
of J. In a sense, J is an “eigenvalue:'’ for a given value of 
radius a (or potential <p 0 ), a solution exists for only one /; this 
J determines the steady corona current I 0 . In general, it is the 
condition (5b) which fixes J. When 7 = 0, equation (9) 
becomes Townsend’s criterion for corona onset; it determines 
the value of a (or <p 0 ), at which the corona becomes 
self-sustained. 

The solutions for 7T_ and ir e which match to equation (7b) 
are 


ic e = (l + y)e~ a a r 
tt_ =(1 + 7)(1 -e~ a a r ) 


( 10 ) 


Equation (9) was used to arrive at these forms. Equation (8a) 
can be rewritten using equation (2) as 


The condition that <p matches to the inner solution (7a) also 
requires that 1 
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Equation (13) with equation (12) determine the value of A as 
a function of 7. After substituting this A (7) into equation (9), 
equation (9) determines 7. 

Equation (12) simplifies if electron attachment is either slow 
(a a = 0) or fast ( a a = oo). In the former case, 


E= + y)r 2 + /4 2 ] 172 


(14) 


and the integral in equation (13) can be evaluated in closed 
form. 

Figure 1 shows A as a function of 7 computed from equations 
(13) and (14). Figure 2 shows a computed from equation (9) 
with the functions A (J) shown in figure 1. One sees in each 
case that a maximum value of a exists above which a corona 
discharge cannot be self-sustained. A maximum 7 exists above 
which the electrode radius becomes zero. One should not be 
mislead by this prediction of a maximum 7: because 
a = a*Ei/<p 0 , when a* is fixed the maximum 7 would occur 
as <p 0 -~ oo. From equation (3b) the corresponding I Q tends to 
infinity like <pl. 


A Strictly a should be treated as a small parameter and A and J should be 
expanded in terms of e; however, the present treatment is accurate to 0(e). 
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Figure 1 .—Coefficient A of equation (7a) as function of current 7(1 +7) for 
7 = 0.01. 



Figure 2.— Scaled electrode radius a versus current 7. For fixed dimensional 
radius a*, "a represents one over the applied voltage. 


Stability of the Corona 

The stability to infinitessimal disturbances of the previous 
solution can be examined by introducing dn$/dt on the left 
side of equations (3a) and forming the linear perturbation 
equations. If the perturbation has the form e st F(r) where 
5 = 0(1), then the asymptotic stability equations can be solved 
in closed form for the case of equation (14). One finds that 
the eigenvalues s all lie in the left-half complex plane, 
indicating stability with respect to this type of exponentially 
growing, axisymmetric disturbances. At present, the stability 
to other types of perturbation is being studied in hopes that 
it will provide insight into corona pulsations. 
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